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f{R) gravity models belong to an important class of modified gravity models where the late time cosmic 
accelerated expansion is considered as the manifestation of the large scale modification of the force of grav- 
ity. f{R) gravity models can be expressed in terms of a scalar degree of freedom by redefinition of models 
' variable. The conformal transformation of the action from Jordan frame to Einstein frame makes the scalar 

degree of freedom more explicit and can be studied conveniently. We have investigated the features of the scalar 
, degree of freedoms and the consequent cosmological implications of the power-law (^J?") and the Starobinsky 

' (disappearing cosmological constant) f{R) gravity models numerically in the Einstein frame. Both the models 

5_( ' show interesting behaviour of their scalar degree of freedom and could produce the accelerated expansion of the 

, Universe in the Einstein frame with the negative equation of state of the scalar field. However the scalar field 

potential for the power-law model is the well behaved function of the field, whereas the potential becomes flat 
for higher value of field in the case of the Starobinsky model. Moreover, the equation of state of the scalar field 
CO , for the power-law model is always negative and less than — 1/3, which coixesponds to the behaviour of the dark 

' energy that produces the accelerated expansion of the Universe. This is not always the case for the Starobinsky 

model. At late times Starobinsky model behaves as cosmological constant A as behaves by power-law model 
^ , for the values of n — > 2 at all times. 
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Although the General Relativity (GR) is a very well-tested and widely accepted theory of gravity, the explanation of the 
late time accelerated expansion of Universe jlll is out of scope within it's ambit. However, because of all the successful test 
performed so far on GR for the low scales, attempts have been made to modify GR in order to account the late time behaviour of 
■ the cosmic expansion. There are two main conceptual approaches that have been led to the modification of GR. In one approach 
00 I a new scalar degree of freedom is incorporated in the energy-momentum tensor on the right hand side of the Einstein's equation, 
, which is dubbed as dark energy for its exotic behaviour with large negative pressure |l2|-Hl- Initially cosmological constant was 
considered as the source of dark energy which was faced with fine tuning problem of an acceptable level. In recent years, as 
, an alternative to the cosmological constant, a variety of scalar field models is proposed to provide a viable explanation for the 
CO ■ phenomenon of late time cosmic acceleration iH-Ht]. In the other approach, the left hand side of the Einstein's equation can be 
modified by considering the late time cosmic accelerated expansion is due to the large scale modification of the force of gravity, 
which leads to the modified gravity theories iflolfTTIl . It needs to be mentioned that there are many other alternative attempts to 
understand the physical mechanism of this late time behavior of the universe (e.g. see ifisll ). 

The simplest class of modified gravity theories is the f{R) gravity theories in which Einstein gravity is modified by replacing 
the Ricci curvature scalar R by an arbitrary curvature function f{R). In the special case when f{R) R, f{R) gravity converge 
to GR without a cosmological constant ifiol [Till . In last few years the f{R) gravity models have been studied extensively in 
the cosmological and astrophysical aspects because of a number of very interesting results (see lITill and see references there 
in). There are many f{R) gravity models which could produce the late time cosmic acceleration, but all of them are not 
cosmologically viable and many of them also suffered from the singularity and stability problem lfl5l[l6|] . So certain restrictions 
have to be imposed on f{R) gravity models to be linearly stable and cosmologically viable. Some of such models can be found 

in iB [li- 
lt is interesting that a new scalar degree of freedom (sometimes referred as scalaron) appears in f{R) theories of gravity due 
to redefinition of models variable. The conformal transformation of the metric to the Einstein frame makes it explicit in the 
action Hol [Til [l6l [isll . So it is possible to extend f{R) gravity to generalized Brans-Dicke (BD) theory with a field potential 
and an arbitrary BD parameter In f{R) gravity the scalar field is coupled to non-relativistic matter (dark matter, baryons) with a 
universal coupling constant — 1/V6, and hence as far as the scalar field is concerned it is convenient work on f{R) gravity in the 
Einstein frame, in which a canonical scalar field is coupled to non-relativistic matter However, conventionally, the Jordan frame, 
where the baryons obey the standard continuity equation {pm oc a^^), is considered as the physical frame in which physical 
quantities are compared with observations and experiments. We are treating essentially the same physics in both frames but with 
different time and length scales that give rise to the apparent difference between the observables in two frames. This apparent 
difference can be overlooked by going back to the Jordan frame when we work on the Einstein frame for some convenience IToll . 

In this work we have studied the ^i?" and the Starobinsky disappearing cosmological constant f{R) gravity |[T2[l models using 
their scalar degree of freedoms in the Einstein frame in the Friedmann-Lemaitre-Robertson- Walker background. The basic aim 
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of this study is to see numerically the features of the scalar field coupled to these f{R) gravity models and hence to see the 
consequent evolution of the cosmological parameters in the Einstein frame. We have considered these two models because: (i) 
The f{R) = ^-R" model has exact power-law solutions and as a result the singularity is not appear in it as suffered by other f{R) 
models. This is manifested in the scalar degree of freedom as the scalar field potential of this model is well behaved function of 
the field ifTgll . (ii) The Starobinsky f{R) gravity model with disappearing cosmological constant fV2i\ is a viable cosmological 
model within the framework of modified gravity which has been studied extensively in last few years in different cosmological 
and astrophysical context |[l4i[l6ll20ll . This will help us to compare the results of ^i?" gravity model with results of the standard 
model IT2!]. We organized this work as follows: In Sec. II we discuss the basics of formalism of the f{R) gravity theory in the 
Jordan frame and then we transform the formalism to the Einstein frame using the conformal transformations to developed field 
equations of the scalar degree of freedom. The numerical analysis of the work based on the formalism developed in the Sec. II is 
done in the Sec. Ill for our models as mentioned above. Finally in the Sec. IV we concludes the results of our analysis. 



II. f{R) GRAVITY FORMALISM IN EINSTEIN FRAME 

In this section we will develop the field equations of f{R) gravity in Einstein frame starting from the action of f{R) gravity 
in the Jordan frame. The f{R) gravity action in the Jordan frame is given by lHoilTTIl . 

(1) 



where = SttG, g is the determinant of the metric g^j^, and Cm is the matter Lagrangian, which is the function of the metric 
g^i, and the matter fields ^'„,. The variation of the action (HJ with respect to g^j^^y the leads to the equation of motion ifiolflilfisll : 

/'-Rm- - V^,/' + (af - i/^ g^, = K^T^,, (2) 

where (/) denotes the derivatives with respect to R. The term □/' does not vanish in modified gravity, which means that there is 
a propagating scalar degree of freedom, defined as, (f — /', whose dynamics is governed by the trace of the equation (|2) given 

by. 



30/' + /'i? - 2/ K^T. 



(3) 



In practice, it is usually difficult to invert the definition of the scalar degree of freedom explicitly for a given f{R) and only 
convenient to determine effective potential in a parametric form lfl6ll . However, using the conformal transformation to the 
Einstein frame it is possible to express the scalar degree of freedom of a given f(R ) in the explicit form fl6ll . In view this we 
now switch to the Einstein frame under the following conformal transformation lUOll : 



_ q2 



(4) 



where il^ is the conformal factor and the quantities in the Einstein frame are represented by a tilde over them. In these two 
frames the Ricci scalars have the following relation: 



R^n^{R + 6Dlu - eg^^'df.Lod^io), 



(5) 



where 



(6) 



For convenient we rewrite the action ([T} in the form |llO| 



5= / d^x^g[—f'R-U 



where 
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Using Eq. ^ and the relation ^/—g = '^V—g, the action O can be written as IT 
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The integral J d'^x^/—gduJ can be made to vanish on account of the Gauss's theorem by using equation (|6]l. Now if make a 
choice ifiotl . 



for /' > and introduce a new scalar field (j) defined by 

K0= ^3/2 In/', 
then the action in the Einstein frame can be expressed as ifioll 



where 



is the potential of the field <p. Here the conformal factor is 

VL^ ^ f ^ exp{y^K(t)), 
which is field dependent. The Lagrangian density of the field ((> is given by 



(10) 



(11) 



(12) 
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From the matter action (fTZt it is clear that the scalar field is directly coupled to matter in the Einstein frame, which can be 
seen more explicitly if we take the variation of the action ([T2|i with respect to the field 4> ifioll - However, we are not interested 
in this context, instead we ignore the matter field to derive the scalar field equations without coupling to the matter Thus taking 
the variation of the action (fT2l i with respect to the field 0, we obtain the field equations without the presence of matter as ifToll 



(16) 



(17) 



Here the dots over (j) indicate the derivatives with respect to conformal cosmic time (i.e. cosmic time in Einstein frame) i. H = ^ 
is the Hubble parameter in the Einstein frame, a{t) being the scale factor in the same frame. If we consider the spatially flat 
Friedmann-Lemaitre-Robertson- Walker (FLRW) spacetime with a comic time-dependent scale factor a{t) and a metric, viz.. 



ds'^ = g^i.dx^'dx'' = -di^ + a^{t)dyi^ 
in the Jordan frame, then the metric in the Einstein frame may be given as 

ds^ = n^ds^ = f'{-df + a^(t)dy.^) = -dP + d'^{t)dx.^, 
which leads to the relations of the cosmic time and scale factor in two frame as 

di = ^/^dt, 5 = y^a, for VT' > 0. 
Finally we define the energy density and the pressure of the scalar field as ([k 



(18) 



(19) 



(20) 



(21) 



(22) 



It is evident from the above formalism that in Einstein frame the scalar degree of freedom for a given f{R) model can be studied 
explicitly and conveniently. In the next section we apply the results of this formalism to study our models of interest. 
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III. SCALAR FIELD DYNAMICS OF / (i?) GRAVITY MODEL 

Using the formalism of the previous section, in this section we will study numerically different features of the scalar field and 
consequent cosmological implications of the power-law and the Starobinsky f{R) gravity models as follows; 



A. Power-law f{R) gravity model 



The general power-law f{R) gravity model is given by 

f{R) = ei?", (23) 

where ^ and n are model parameters. We use this f{R) gravity model in this study because of its unique character for the 
existence of power-law solutions and hence absence of the singularity problem as mentioned in the Sec. I llT9ll . Here our main 
emphasis is to see the expUcit behavior of scalar degree of freedom possessed by f{R) gravity model and its usefulness to study 
the cosmological evolution. 

Using the equation ( fTTT i. the scalar field (j) corresponding to the model (|23] | can be expressed as 

ll^-r ^ ^-llni^nR--^). (24) 
Using this expression for (f), the field potential ( fTSi t takes the form 

2-n 

^ [ ^^ ) ■ ^'^^ 

This equation for the field potential shows that the value of n should be within two possible ranges: (i) ri < 1 and (ii) 1 < n < 2 
for a viable field potential. For our numerical simulation in this work we consider the second range of the values of n and the 
values of ^ < 1 . We consider k ~ 1 through out this work. 

Fig. [T]shows the variation of field potential with respect to the field (f>. The left panel is for n = 1.5 with ^ — 0.1, 0.3 and 
0.5, whereas the right panel is the plot for n = 1.5 and 1.6 with ^ = 0.1. From this figure it is seen that the potential is well 
behaved function of the field and the potential is small for negative value of the field and it increases very fast as soon as the 
field (f) become positive. This process is slowing down when the value of ^ increases, which is more effective for lower value 
range of ^, as clear from the left panel. On the other hand as seen from the right panel that there is no significant difference of 
the potentials for the different values of n within the range of negative field, but the difference becomes more prominent as soon 
as the field (f> acquired the positive values. The effect of n on the value of the potential is similar to that for ^, but the effect on 
higher positive value of (f> is more noticeable in this case than the case for ^. 




FIG. 1: The variation of field potential with respect to the field 0. The left panel is for n = 1.5 with ^ — 0.1, 0.3 and 0.5, whereas the right 
panel is for n — 1.5 and 1.6 with ^ = 0.1. The value of k is taken as one. 

Using the potential (l25l l. the coupled field equations ( fT6] l and (fTTI i are solved numerically for different values of n and ^ as 
mentioned above. The results of the numerical solution are shown in the Fig. |2] The left panel shows the numerical solutions 
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for n = 1.5 with ^ = 0.1, 0.2 and 0.3. The middle panel is for n = 1.5 and 1.6 with ^ = 0.1. In both cases the initial field 
value and the initial field velocity are taken as 1.0 and 0.1 respectively. As clear from the left panel that the value of the field (f> 
falls towards to the negative values starting from its initial value as time increases. For the initial period of cosmic time cf) falls 
very fast and this tendency slows down as time passage. For higher values of ^ the overall tendency of the to fall towards the 
negative values is less in comparison to its lower values. This is similar to the nature of the field potential with respect to ^. 
Form the middle panel it is observed that the field falls faster initially for the smaller value of n as compared to the higher value, 
but the situation reversed in the latter cosmic times. During the time before the situation is reversed, the difference of the field (j) 
is not very high, however after this period the difference become increasingly significant as time elapsed. We have also studied 
the effect of the initial field and field velocity on the time evolution of the field, whose results are shown in the right panel of the 
Fig. 121 This plot is for n = 1.5 and ^ = 0.1 with (po = 1.0, 0.1 and (po = 0.1, 0.01, 0.001. We have seen that the initial values 
of <j> and (f> does not have much noticeable effect on the time evolution of 0. Initial value of cj) has slight effect during initial 
period, which eliminates gradually in latter times. On the other hand the initial values of (/> does not have any effect on the time 
evolution of cj). 




FIG. 2: Numerical solutions of equations J16t and ( I17t using the field potential l l25i l for different values of n and ^. The left panel is for 
n = 1.5 with different values of ^. The middle panel is for ^ = 0.1 with two values of n. In both these panels <j)o — 1.0 and (po = 0.1 are 
used. While the right panel is for different initial values of and with n = 1.5 and ^ = 0.1. The value of k is taken as one. 



From these numerical solutions we have calculated the scale factors d{t) in Einstein frame for different values of n and ^ as 
mentioned above. The time (cosmic) evolution of the scale factors obtained from these calculation are shown in the Fig. |3] The 
left panel of this figure shows the time evolution of the scale factor for n = 1.5 with different values of ^. Whereas right one is 
for ^ = 0.1 with two different values of n. The initial parameters are same as for the first two panels of the Fig. |2] It is observed 
from this figure that the f{R) = ^i?" gravity model leads the very fast expansion of the Universe in the Einstein frame (will also 
lead similar expansion in the Jordan frame) similar to the exponential expansion. The rate of expansion of the Universe slowed 
down gradually for higher values of ^ in comparison to its lower values. During initial period this process is slow, but become 
more significant as time passage. On the other hand the expansion is slower for lower value of n in comparison to the case for 
the higher value of n as clear from the right panel. In the initial period the difference of the values of the scale factor is very lass, 
but becomes more prominent in the latter times. 

At last we have calculated the equation of state w ~ jP^ in Einstein frame using equations (ISTT i and (l22T i from the numerical 
solutions mentioned above, which is shown in the Fig. [3] Panels of this figure are for the same model parameters respectively 
as for the panels of the Fig. [3] It is observed from the figure that the equation of state w of the scalar field <p for the model ( |23] l 
always less than —0.75 for all cases of our study and for all times. For higher values of ^, w is smaller than that for smaller value 
of ^ during period for some initial times, but for latter times w becomes equal for all ^. But for higher value of n, w is always 
smaller than that for the smaller value of n with a similar trend with respect to time and ^ as mentioned above. As the equation 
of state is always negative, therefore the model (|23] l could produced the accelerated expansion of the Universe and it corresponds 
to dark energy, whose equation of state is w < —1/3 (in Jordan frame, by conformal transformation from Einstein frame to the 
Jordan frame we can show that in our model also the equation of state lies within this range). If we increase the values of n then 
the equation of state w in Einstein frame will approach to —1. This implies that the model (l23T l with higher values of n (nearly 
equal to 2) corresponds to the cosmological constant A. 
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FIG. 3: Evolution of the scale factor in the Einstein frame obtained from the numerical solutions of equations ( 116) and ( |17t using the field 
potential i25\ for different values of n and ^. The left panel is for n = 1.5 with different values of ^. Whereas the right panel is for ^ — 0.1 
with two vales of n. In both these panel (t>o = 1.0 and (j>o = 0.1 are used. The value of k is taken as one. 




FIG. 4: Equation of state of the scalar field in the Einstein frame obtained from the numerical solutions of equations ( 116) and l ll7t using the 
field potential i25\ for different values of n and ^. As in the previous figure the left panel is for n = 1.5 with different values of ^. Whereas 
the right panel is for ^ = 0.1 with two vales of n. All initial conditions and model parameters are same with the previous figure. 



B. Starobinsky f{R) gravity model 

The Starobinsky ,f{R) gravity model with disappearing cosmological constant lfl2ll is. 



fiR) = R + XRQ 



R^ 



- 1 



(26) 



where n, X and Rq are free parameters. The values of n and A are greater than zero, whereas the value of Rq is of the order of 
the presently observed cosmological constant A(= STrGpvac)- This is carefully constructed model to give viable cosmology and 
to be satisfied with solar system and laboratory tests. Notwithstanding it suffers from the singularity problem lfl6l . which could 
be cured by adding by a cx i?^ term to it fl^ . 

The scalar field embodied in this model (l26T l which is obtained by using the equation ( fTTT ) can be written as 



31 
22 



Inf 



3 1 

22^ 



In 



1 - 2Xn 



R 
Ro 



-{"+!)■ 

Rl 



(27) 



7 



If we use this expression for (/> in a region where R/Rq >> 1, the field potential (flJl i for this model becomes, 



(2n + 1) 



1 — exp{yj'2/3K 
2nA 



(28) 



This expression for the field potential shows that for the validity of the Starobinsky model Goi the values of n and A must 
be greater than zero, but may take any values freely beyond it in contrast with the power-law model (l23T l. However, for our 
numerical work on this model, we have taken n = 1 and 2 with the values of A = 1.2, 1.5, 2.0 and 4.0. We have followed the 
same numerical procedures for this model also that have been used for the case of power-law model. 

The behaviour of the scalar field potential (|28^ with respect to the field (j> associated with the Starobinsky model (l26T l for 
different values of n and A as mentioned above is shown in the Fig. |5] We observed that for smaller values of A, when the field cj) 
is rolling from negative to zero, the potential increases very fast from negative to its maximum positive value as the field reached 
its zero value and then falls off rapidly to a almost flat region near its zero value when the field increases again from its zero 
value. This behaviour of the potential for such A values form a pattern like the A-pattem. Beyond this range of smaller values 
of A, the A-pattem of the potential vanishes gradually and the potential always remain positive for all values of the field (p. The 
range of this smaller values of A depends on the value of n, higher the value of n, smaller is this range. Thus for higher values of 
n and A, the potential falls off from a very high value to its almost flat region near zero. Moreover, for all values of ji and A the 
potential becomes flat at a particular high positive value of (j> depending on the value of A, higher the value of A, slightly higher 
this value of ch. 




^ ^ <t> 

FIG. 5: The behaviour of the field potential ( I28t with respect to the field (j> corresponding to the Starobinsky model i26\ . The left panel is for 
n — 1 and the middle panel is forn = 2 with A = 1.2, 1.3, 2.0 and 4.0. The right panel is for n — 1 and 2 with A — 1.5 to show exclusively 
the effect of the parameter n on the potential for a given value of A. 



The Fig. |6] shows the results of numerical solution of field equations ( fTSI l and ( fTTI i using the field potential ( [28] ) for different 
values of n and A. The initial values of field <j> and field velocity (f> are same as their values used in the case of power-law model. 
The general pattern of time evolution of the field is that the field cj) falls off rapidly from its initial positive value to its minimum 
negative value and then oscillates for sometimes with a decreasing amplitude before becomes steady with a negative value. For 
higher values of A and n the field is legging behind the field for the lower values of these parameters by shifting its minimum 
towards more negative side. This indicates that the magnitude of the field is smaller for higher values of A and n in comparison 
to the field of their smaller values. It should be noted that for the given A or the given n this behaviour of the field (f> is observed 
for different values of n or A. However, in the case of A = 4.0 for n = 2 a slight different behaviour of the field is observed 
as seen from the top right panel of this figure wherein the field falls from higher initial values in contrast with the behaviour for 
other cases. The initial field velocity has noticeable effect on the time evolution of the field as it is observed from the bottom 
right panel of this figure that for the lower value of the initial field velocity the field magnitude increases for a period before the 
field (f) becomes steady with the same value attained for other initial conditions. 

The study of the scale factor evolution in the Einstein frame using the scalar degree of freedom of this model as we have studied 
in the case of the power-law f{R) model is shown in the Fig. |7] It is clear from this figure that during a considerable initial 
period the expansion of the Universe is slower and after that period the Universe experiences very fast accelerated expansion 
similar to the exponential expansion under this model. This initial period becomes shorter as the value of A increases. Moreover 
for a given value of A this slower expansion period of the Universe is shorter for the higher value of n as shown in the right panel 
of this figure. That is higher the value of n and A the earUer is the period of accelerated expansion of the Universe. Thus this 
model could explain the late time accelerated expansion of the Universe with the suitable values of n and A. 
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t t 

FIG. 6: Time evolution of the field (j> obtained from numerical solutions of equations il6\ and ( I17t using the field potential i28t for different 
values of n and A. The top two panels are for n = 1 (left) and 2 (right) with different values of A. The bottom left panel is for A = 1.2 with 
two different values of n. In all these panels (po = 1.0 and 00 = 0.1 are used. While the bottom right panel is for different initial values of (j) 
and <j) with n = 2 and A = 2.0. 

As described in the previous Sec.III(A) for the power-law / {R) gravity model we have calculated in the same way the equation 
of state of the scalar field related with the Staroninsky f{R) gravity model ( [26] l. The results of this numerical calculation is shown 
in the Fig. [8] The patterns of time evolution of the equation state of the scalar field for different values of A and n are almost 
similar to the patterns of time evolution of the field for theses parameters. That is the equation of state is also oscillate with 
decreasing amplitude during a period after falling from its initial positive value and then after that period it remains steady with 
time. During a initial period depending on the value of A and n the equation of state falls from 1 to —1, which indicates that 
during this period the scalar field corresponding to the model (|26] | make a transition from free scalar field {w = 1) through 
matter (w « 0, radiation (w = 1/3) and dark energy (w < —1/3) to cosmological constant A (w = —1). During the oscillation 
period the scalar field behaves as matter to the cosmological constant and finally behaves as the cosmological constant A after 
completing this period depending on the values of the parameters A and n. There is a special behaviour of the equation of state 
for the n = 2 and A = 4.0 during the initial period. In contrast to the behaviour of the equation of states for the other pairs of 
parameters, the equation of state of this pair of parameters initially make transit from the dark energy state to other states and 
then make transition to the cosmological constant state. However for latter periods the behaviour of equation of state for this 
pair of parameters is similar to the other cases. 

IV. CONCLUSION 

We have investigated the scalar degrees of freedom and it's related cosmological implications of two f{R) gravity models, 
viz., (i) the power-law f{R) gravity model ^i?" and (ii) the Starobinsky disappearing cosmological constant f{R) gravity 
model for different values of the parameters of these models. This study is basically on the theoretical interests rather than the 
observational consequences. The power-law f{R) gravity model is considered due to the existence it's power-law solutions and 
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FIG. 7: Evolution of the scale factor in the Einstein frame obtained from the numerical solutions of equations ( 116) and ( |17t using the field 
potential ( 128b for different values of n and A. The left panel is for n = 1 with different values of A. Whereas the right panel is for A = 2.0 
with two vales of n. In both these panel (f)o = 1-0 and (j>o = 0.1 are used. 




t t t 

FIG. 8: Equation of state of the scalar field in the Einstein frame obtained from the numerical solutions of equations ( 116) and l ll7t using the 
field potential ( 128b for different values of n and A. The first two panel is for n = 1 and 2 respectively with different values of A. Whereas the 
last panel is for A — 4.0 with two vales of n. All initial conditions are same as in the previous figure. 



consequently absence of singularity problem liT9ll . On the other hand the Starobinsky f{R) gravity model with the disappearing 
cosmological constant lfl2ll is a standard f{R) gravity model which provide viable cosmology by satisfying the solar system and 
the laboratory tests. As the scalar degree of freedom of f{R) gravity model can be treated explicitly in the Einstein frame, so 
we used this frame for this study by conformal transformations from the Jordan frame. 

It is found that the field potential of the power-law f{R) gravity model is the well behaved function of the the field (j) 
representing the scalar degree of freedom of the model and increases very fast as 4> increases, which depends on the model 
parameters ^ and n such that for the lower value of them such process is faster than higher value. From the numerical simulation 
of the field equations, we found that (p falls towards the negative values as time increases, which is more faster during the initial 
period of time than late times. The rate of falling the value of the field slower for the higher values of the parameter ^. However 
n does not have much impact as ^ has on the behaviour of field cj) with respect to time. Initial conditions on the field and field 
velocity have negligible effect on the nature of the time evolution of (f>. 

In the case of the Starobinsky f{R) gravity model with disappearing cosmological constant the field potential shows a peculiar 
behaviour with respect to the field (j). When the field increases from its negative value to zero the field potential increases from 
its lowest negative value to the maximum positive value and then falls off very fast to the flat region near to its zero value as the 
field increases from zero to higher positive values for a lower value range of the parameter A. This range of the lower value of the 
parameter A decreases with increasing the value of n. This behaviour of the filed potential form a pattern similar to the A-pattern. 
The A-pattern gradually vanishes with increasing values of the parameters A and n. For higher values of the field cj) the values 
of field potential for different A and n parameters are almost identical with a slight dependence on the value of the parameter A. 
The numerical solution for the time evolution of the field for this model shows that the field falls off rapidly from its maximum 
positive value to its minimum negative value in the earlier period and then starts oscillating with a decreasing amplitude with 
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time for a period before become steady with time taking a particular negative value. This pattern of time evolution of the field 4> 
depends on the the model parameters in such a way that with increasing values of A and n the field (f) leg behind its values for 
lower model parameters with decreasing amplitudes to attain steady value much earlier 

The study of the time evolution of the scale factor in the Einstein frame shows that the Universe is accelerating very fast for 
both of our f{R) gravity models, which depends on values of the parameters of the models. For the case of the power-law f{R) 
model the scale factor depends on the parameters ^ and n in such a way that accelerating process slows down for higher values 
of ^ and for lower values of n. Whereas in the case of the Starobinsky f{R) model the accelerating process becomes more faster 
for the higher value of both parameters A and n. Moreover, in this case there is a initial period depending on the value of A in 
which the expansion of the Universe is very slow. This initial period is shorter for the higher value of the parameter A and does 
not depends on the parameter n. After this initial period the Universe is expanding with a very fast acceleration process. This 
indicates that this model with the suitable model parameters A and n is very effective to give the late time accelerated expansion 
of the Universe. 

The equation of state of the scalar field for the power-law model is always negative and less than —1/3, which corresponds to 
the behaviour of the dark energy that produces the accelerated expansion of the Universe. The equation of state of the field for 
this model tends to —1 as the value of n increases towards higher values, nearly equal to 2. Thus for such values of n, the scalar 
field for the model behave as the cosmological constant A. The same study for the scalar field of the Starobinsky model shows 
that the equation of state oscillates with a decreasing amplitude with time for a certain period after falling from its maximum 
value around 1 to its minimum value at —1. After this period the equation of state takes the value equal to —1. This implies that 
the scalar field related with the Starobinsky disappearing cosmological constant f{R) gravity model passages through different 
phases before finally behaves as the cosmological constant A. The behavior of the time evolution and the dependence of the 
equation of state on the model parameters A and n are almost similar to the field itself. 
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